Monte Carlo reconstruction of the inflationary potential 
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We present Monte Carlo reconstruction, a new method for "inverting" observational data to 
constrain the form of the scalar field potential responsible for inflation. This stochastic technique is 
based on the flow equation formalism and has distinct advantages over reconstruction methods based 
on a Taylor expansion of the potential. The primary ansatz required for Monte Carlo reconstruction 
is simply that inflation is driven by a single scalar field. We also require a very mild slow roll 
constraint, which can be made arbitrarily weak since Monte Carlo reconstruction is implemented 
at arbitrary order in the slow roll expansion. While our method cannot evade fundamental limits 
on the accuracy of reconstruction, it can be simply and consistently applied to poor data sets, 
and it takes advantage of the attractor properties of single-field inflation models to constrain the 
potential outside the small region directly probed by observations. We show examples of Monte 
Carlo reconstruction for data sets similar to that expected from the Planck satellite, and for a 
hypothetical measurement with a factor of five better parameter discrimination than Planck. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Given a detailed model of physics at the Grand Unified 
Theory (GUT) scale that supports inflation one can de- 
duce many of the overall properties of the universe from 
first principles. However, new strides in precision cos- 
mology are likely to present us with the inverse problem: 
given an accurate measurement of the primordial pertur- 
bation spectrum, can we determine the particle physics 
that drove inflation? This question leads to "reconstruc- 
tion" - rebuilding a portion of the inflationary potential 
from observations - whose apotheosis is the work of Lid- 
sey et al. Q . The conclusion of this effort was that while 
a limited reconstruction may be possible, it is overly op- 
timistic to hope to deduce more than the rough form 
of a small piece of the inflaton potential on the basis of 
astrophysical observations alone. 

Despite this apparent setback, the prospect of being 
able to probe ultra-high energy particle physics via as- 
trophysical measurements remains tempting. In this pa- 
per we present a new approach to the task, Monte Carlo 
reconstruction. Unlike previous workers, we tackle the 
problem stochastically, by generating large numbers of 
inflationary models and their associated potentials and 
identifying those that lead to a spectrum of primordial 
perturbations which fits inside a specified window of pa- 
rameter space. We sacrifice the prospect of unambigu- 
ously reconstructing the inflationary potential, but in- 
stead map out the ensemble of viable potentials. We can 
then assess the extent to which these potentials resem- 
ble one another, and thus the degree of ambiguity in the 
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"reconstruction" . 

The starting point for Monte Carlo reconstruction is 
the Hubble slow roll expansion 0, which leads to a set 
of first order "flow" equations ]3|, Q that describes the in- 
flationary dynamics. The Hubble slow roll expansion can 
be systematically continued to arbitrary order but, as we 
explain below, the truncated expansion actually yields an 
exact solution of the full equations, drawn from a subset 
of the overall solution space. We select a "universe" by 
assigning random values to each term in the truncated 
slow roll expansion, and integrate the flow equations for- 
ward to the end of inflation (or far enough into the fu- 
ture to show that the solution approaches an inflationary 
attractor). We can then determine the epoch that corre- 
sponds to cosmological structure formation and thus the 
predicted parameters of the primordial spectrum, allow- 
ing us to select the models which fit into a pre-determincd 
window of parameter space. 

As with any other reconstruction program, we face lim- 
itations dictated by the physics of inflation itself. Since 
the field rolls slowly, it traverses only a small portion of 
the potential during the epoch in which the primordial 
perturbations are laid down. Consequently, even if all the 
potentials we recover for a specific set of cosmological pa- 
rameters overlap in the region where the cosmologically 
relevant portion of the spectrum is laid down, they will 
still diverge outside this region. However, because we 
reconstruct the entire potential, we can estimate the de- 
gree to which the reconstructed potentials diverge from 
one another. 

In the following section we summarize the inflationary 
dynamics and the flow equations, and in Section 3 ap- 
ply these to spell out the process of Monte Carlo recon- 
struction in detail. We present the results of a variety of 
simulations in Section 4, and test how accurately Monte 
Carlo reconstruction can recover a known inflationary po- 
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tential. As expected, a fully consistent reconstruction 
requires the unambiguous detection of tensor modes in 
the CMB, but even in the absence of a detectable tensor 
signal we can still recover some information about the 
inflationary potential. If we assume the error bars on the 
scalar spectral index and tensor-to-scalar ratio that are 
expected for the Planck mission, the underlying potential 
cannot be unambiguously recovered. However, if we de- 
crease these error bars (perhaps by adding information 
from large scale structure surveys, or better measure- 
ments of CMB polarization) by a factor of 5, then we 
can start to recover a 4> n potential. Our conclusions are 
summarized in Section 5. 



II. INFLATION AND FLOW 

We consider inflation driven by a single homogeneous 
scalar field <fr (the inflaton) with potential V(</>) and equa- 
tion of motion 



4> + 3H^ + V'(^) = 0, 



(1) 



where H = (a/ a) is the Hubble parameter. We assume a 
spatially flat, homogeneous and isotropic universe where 
the inflaton field is the only contribution to the energy- 
momentum tensor, so that the Einstein field equations 
have the familiar form 



H z 



377Jpj 



v(4>) + \ft 



and 



8tt 
3mpj 



vw - if 1 



(2) 



(3) 



Here m P1 = G -1 / 2 ~ 10 19 GeV is the Planck mass. These 
background equations, along with the equation of motion 
(|l|) , form a coupled set of differential equations that de- 
scribe the evolution of the universe. The limit = 
corresponds to a de Sitter universe, with the scale factor 
increasing exponentially in time 
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const, 



(4) 



In general H is not exactly constant, but varies as the 
field <fi evolves along the potential V(4>). A powerful way 
of describing the dynamics of an inflationary universe 
with a varying field (and non-trivial potential) is to ex- 
press the Hubble parameter as a function of the field </>, 
H = H(<f>), which is consistent provided <f> is monotonic 
in time. The equations of motion become Jl], ||, ||, (jj 
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These are completely equivalent to the second-order 
equation of motion (Q). The second of the above equa- 
tions is referred to as the Hamilton- J acobi equation, and 
can be written in the useful form 
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where e is defined to be 
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The physical meaning of e can be seen by expressing Eq. 
(0) as 
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(8) 



so that the condition for inflation (a/a) > is given by 
e < 1. The scale factor is given by 
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where the number of e-folds ./V is 
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It is convenient to use N as the measure of time during 
inflation. We take t e and <\> e to be the time and field 
value at end of inflation. Therefore N is defined as the 
number of e-folds before the end of inflation and increases 
as one goes backward in time, dt > dN < 0. The 
sign convention for y/e must be applied carefully, and we 
take it to have the same sign as H'(<fr): 
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Liddle, Parsons and Barrow use this as the starting point 
for an infinite hierarchy of slow roll parameters : 
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The evolution of these parameters during inflation is de- 
termined by a set of "flow" equations pi, 0] , 
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Hero the time variable is the number of e-folds N, where 
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(14) 



The derivative of a slow roll parameter at a given order is 
higher order in slow roll. Taken to infinite order, this set 
of equations completely specifies the cosmological evolu- 
tion, up to the normalization of the Hubble parameter 
H. When truncated at finite order, by assuming that 
the ^Ah are all zero above some fixed value of £, the solu- 
tion of the flow equations still yields an exact solution of 
the background equations, albeit one that is drawn from 
a special subset of the overall solution space. 



III. MONTE CARLO RECONSTRUCTION 

The inflationary dynamics are fully specified by the 
values of the slow roll parameters [e, a, Ah] at a fixed 
time, which serve as initial conditions for the flow equa- 
tions, truncated at order M in the slow roll expansion. 
We now show that this information determines the infla- 
tionary potential V(<p), up to a constant multiplier. The 
starting point is the Hamilton-Jacobi equation, 
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We have e(N) trivially from the flow equations. In order 
to calculate the potential, we need to determine H(N) 
and <fi(N). With e known, H(N) can be determined by 
inverting the definition of e, Eq. (Q): 



1 dH 
HdN 



(16) 



This can be viewed as the lowest-order member of the 
system of flow equations (^). Similarly, (j>(N) follows 
from the first Hamilton-Jacobi equation (0): 
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Using these equations and Eq. (p.5[), the form of the po- 
tential can then be fully reconstructed from the numeri- 
cal solution for e(TV). 1 The only necessary observational 
input is the normalization of the Hubble parameter H, 
which enters the above equations as an integration con- 
stant. Here we use the simple condition that the den- 
sity fluctuation amplitude (as determined by a first-order 
slow roll expression) be of order 10 , 
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1 Expressing the potential in this way is certainly not new 
is the basis of "functional reconstruction" method of Refs. 
and well as the "Stewart-Lyth inverse problem" of Rcf. |l0| 




A more sophisticated treatment would perform a full nor- 
malization to the COBE CMB data gl| |§ . The value of 
the field, <fi, also contains an arbitrary, additive constant. 

Given a solution to the flow equations it is straightfor- 
ward to determine the observational predictions for that 
model, even before the potential is computed. For in- 
stance, measurements of CMB anisotropies |13, [l4| may 
determine the tensor/scalar ratio r, the spectral index n, 
and the "running" of the spectral index dn/dlogk, and 
we focus on these quantities here. To lowest order, the 
relationship between the slow roll parameters and the 
observables is especially simple: r — e, n — I = a, and 
dn/dlogk = 0. To second order in slow roll, the observ- 
ables are given by ||, [lfj , 



r = e[l-C(<r + 2e)], (19) 
for the tensor/scalar ratio, and 

n - 1 = a - (5 - 3C*) e 2 - i (3 - 5C) <re+ 1 (3 - C) ( 2 A H ) 

(20) 

for the spectral index. The constant C = 4(ln2 + 7), 
where 7 ~ 0.577 is Euler's constant. Derivatives with 
respect to wavenumber k can be expressed in terms of 
derivatives with respect to N as plf 
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The scale dependence of n is then given by the simple 
expression 
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which can be evaluated to third order in slow roll by 
using Eq. ( p0| ) and the flow equations. The final result 
following the evaluation of a particular path in M dimen- 
sional "slow roll space" is a point in "observable param- 
eter space", i.e. (r, n, dn/dlogk), corresponding to the 
observational prediction for that particular model. This 
process can be repeated for a large number of models, 
and used to study the attractor behavior of the inflation- 
ary dynamics. In fact, the models cluster strongly in the 
observable parameter space || . 

We are now in a position to ask the key question: given 
an observational constraint, what dynamics - and under- 
lying inflationary potential - are compatible with that 
constraint? We proceed by generating a large number of 
"universes" with random values (in a sense that is made 
explicit below) for the slow roll parameters. Given an 
allowed region in the observable parameter space, cen- 
tered on specified values of r,n, and dn/dlogk with a 
width given by the error bars on these quantities, we 
can construct an ensemble of inflationary potentials con- 
sistent with with the specified "window" in parameter 
space. This Monte Carlo approach to reconstruction dif- 
fers from previously used techniques pi [l^, |l7) in that 
it generates an ensemble of potentials consistent with a 
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given constraint, rather than attempting to produce a 
parameterized form for the potential and then constrain 
those parameters. The advantage of Monte Carlo recon- 
struction is that by evaluating the flow equations to high 
enough order we need only very mild a priori assump- 
tions about the form of the potential. In addition, by 
probing the parameter space in a uniform way, we obtain 
what is in some sense a "fair sample" of the potentials 
consistent with a given observational bound. Even so, we 
do not have a metric on the space of initial conditions, 
so it is not possible to derive statistical inferences from 
the ensemble of potentials. 

The condition for the end of inflation is that e = 1. 
Integrating the flow equations forward in time will yield 
two possible outcomes. One possibility is that the con- 
dition e = 1 may be satisfied for some finite value of N, 
which defines the end of inflation. We identify this point 
as N = so that the primordial fluctuations are actually 
generated when N ~ 50. Alternatively, the solution can 
evolve toward an inflationary attractor with r = and 
n > 1, in which case inflation never stops. 2 In reality, 
inflation must stop at some point, presumably via some 
sort of phase transition, such as the "hybrid" inflation 
mechanism Jl8|, [l|, |2(J. Here we make the simplifying 
assumption that the observables for such models are the 
values at the late-time attractor. 

To summarize, the algorithm for Monte Carlo recon- 
struction is as follows: 

1. Specify a "window" of parameter space: e.g. cen- 
tral values for n — 1, r or dn/dlnk and their asso- 
ciated error bars. 

2. Select a random point in slow roll space, [e, rj, i \n}, 
truncated at order M in the slow roll expansion. 

3. Evolve forward in time (dN < 0) until either (a) 
inflation ends (e > 1), or (b) the evolution reaches 
a late-time fixed point (e = e \n = 0, a = const.) 

4. If the evolution reaches a late-time fixed point, cal- 
culate the observables r, n — 1 , and dn/d\nk at this 
point. 

5. If inflation ends, evaluate the flow equations back- 
ward N e-folds from the end of inflation. Calculate 
the observable parameters at this point. 

6. If the observable parameters lie within the speci- 
fied window of parameter space, compute the po- 
tential and add this model to the ensemble of "re- 
constructed" potentials. 

7. Repeat steps 2 through 6 until the desired number 
of models have been found. 



In principle it is possible to carry out Monte Carlo re- 
construction with no assumptions about the convergence 
of the hierarchy of slow roll parameters. In practice, the 
flow equations (|T^ ) must be truncated at some finite or- 
der and evaluated numerically. Moreover, for any given 
path in the parameter space, we do not know a priori the 
correct number of e-folds N at which to evaluate the ob- 
servables, since this depends on details such as the energy 
density during inflation and the reheat temperature [Q. 
Consequently, after truncating to order M in slow roll, 
we select our models' parameters randomly from the fol- 
lowing uniform distributions: 

N = [40,70] 

e = [0,0.8] 

a = [-0.5,0.5] 

2 A H = [-0.05,0.05] 

3 Ah = [-0.025,0.025], 

M+1 A H = 0. (23) 

and so forth, reducing the width of the range by factor of 
five for each higher order in slow roll. The series is closed 
to order M by taking m+1 Ah = 0. We use M = 5 for the 
calculations in this paper. The exact choice of ranges for 
the initial parameters does not have a large influence on 
the result of the Monte Carlo process, as long as they are 
chosen such that the slow roll hierarchy is convergent. As 
noted above, the form of the flow equations ([D]) ensures 
that the derivative of ^Ah depends only on parameters 
of order t and £ + 1 , so this truncation still leads to an 
exact evaluation of the flow equations to infinite order. 
We are selecting a finite subset out of an infinite number 
of possibilities for initial conditions, but the background 
evolution for a given model is evaluated exactly. 3 

The result of integrating the flow equations for a par- 
ticular model is a "path" in the slow roll space param- 
eterized by the number of e-folds N. Figures 1-3 show 
examples of paths plotted in the a — e plane for different 
assumptions about the central values for the parameters 
r, n, and dn/dlogk. Also plotted is the corresponding 
reconstructed potential for each path. Surprisingly, while 
quite complex behavior is possible for the slow roll pa- 
rameters, this behavior is only weakly reflected in the 
shape of the potential itself. 



IV. MONTE CARLO RECONSTRUCTION IN 
PRACTICE 

We now describe two concrete applications of this 
method. First, we select central values in various re- 
gions of the observable parameter space with the error 



2 See Ref. for a detailed discussion of the fixed-point structure 
of the slow roll space. 



As usual in all analyses of inflation we ignore the back reaction 
of quantum fluctuations on the background evolution. 
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bars expected from Planck, and reconstruct the infla- 
tionary potential based on this "synthetic" data. We 
find that Planck will allow us to determine the quali- 
tative form of the potential, but is insufficiently precise 
for well-constrained reconstruction. This is consistent 
with the conclusions of other studies |Q, [Tq| . Second, we 
choose a region of parameter space centered on the values 
for n, r and dn/dlnk for a specified potential - in this 
case V((j>) oc c^ 4 - and apply our reconstruction algorithm 
with different sized error bars, thus determining the ob- 
servational precision that is needed in order to reliably 
reconstruct a known potential. 

Figure 4 shows the results of two different "reconstruc- 
tions": for r = 0.02, n = 0.95 and dn/d\ogk — 0, and 
for r = 0, n = .93 and dn/dlogk = 0. In both cases the 
errors bars are those anticipated for the Planck mission's 
measurement of the spectrum [^4| p2]| , namely Sr ~ 0.01, 
Sn ~ 0.01 and Sdn/dlogk ~ 0.01. The first case has 
a tensor contribution that would be resolved by Planck, 
and the normalization of the potential is relatively tightly 
constrained, V ~ 1 0~ n mpj. The second reconstruction 
assumes that the tensor amplitude is not detectable by 
Planck, with a central value r ~ 0.0. In this case, the 
normalization of the potential is very poorly constrained, 
although the shape of the potential is consistent with 



standard "small- field" models |L3]]. The clear conclusion 
from these plots, and the other cases that we have exam- 
ined, is that Planck may be able to determine the quali- 
tative features of the potential, but it will not permit the 
quantitative reconstruction of the inflationary potential 
on its own. 

It should not come as any particular surprise that 
Monte Carlo reconstruction does not yield a tight con- 
straint on the possible form of the inflationary poten- 
tial. However, the question this immediately raises is how 
much more accurately - relative to the precision expected 
from Planck - would we need to measure the spectrum in 
order to be able to make a quantitative statement about 
the functional form of the potential. We tackle this issue 
by choosing values of n, r and dn/dlogk consistent with 
a known potential, in this case the canonical A</> 4 model, 
which is typical of "large-field" models and has a tensor 
fluctuation amplitude observable by Planck. With this 
assumption, the central values for the spectral index and 
tensor-to-scalar ratio are n = 0.943 and r = 0.02 (evalu- 
ated 50 e-folding before the end of inflation). For V oc 4 , 
dn/dlogk is close to zero, and we take the central value 
to be precisely zero. 

Using the above choices for the central values, we 
have performed Monte Carlo reconstruction with Planck- 
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sized errors bars and with Sn = 0.002, Sr = 0.002 and 
Sdn/d\ogk — 0.005. The second choice leads to a box in 
the (n, r) plane that is 25 times smaller than the bound 
predicted for Planck and thus represents a substantial 
improvement in experimental precision. We divided the 
dn/dlogk bound by 2 rather than 5 solely for computa- 
tional convenience. The second reconstruction involved 
computing approximately 90 million "trial" inflationary 
models in order to filter out 100 solutions that passed 
through the specified window in parameter space. This 
took several days of CPU time on a fast desktop machine, 
and would have taken even longer if we had applied the 
same scaling to Sdn/dlogk as we did to Sn and Sr. 

If we plot graphs to analogous to those in Fig. 4 for 
these two reconstructions we find that, unsurprisingly, 
the overall shape of both sets of reconstructed potentials 
are consistent, but that there is less spread in the set with 
the smaller error bars. However, to determine whether 
any sort of quantitative reconstruction is possible, we 
adopt the prior that the potentials are proportional to 
(<t> — 4>o) m , where 0o is an unknown offset, and then per- 
form a least-squares fit to determine the best-fit value of 
m. In practice, it is sufficient to normalize the height of 
the reconstructed potential to precisely unity at 4> — 



and then fit it to the following functional form: 

(1 - c4) m (24) 

This choice forces the fitted (and rescaled) potential to 
be unity when = 0. Adding normalization as a third 
parameter to the fit does not have a significant impact 
on the computed values of m. 

The results for both simulations are shown in Fig. 5. 
The error bars expected from the Planck mission do not 
allow one to conclude that V oc tp 4 , but the tighter 
bounds on n and r do rule out values of m markedly 
different from 4. 4 

We thus tentatively conclude that while Planck cannot 
measure the perturbation spectrum accurately enough to 
put even mild quantitative constraints on the potential, 
increasing the accuracy with which both n and r can be 



4 In reality, if the computed value of m is significantly different 
from 4, what we actually learn is that the potential cannot be 
well described by Eq. - since setting m to a value signifi- 

cantly different from 4 produces values of n and r well outside 
our assumed range. 
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FIG. 3: Paths in the (a, e) plane (left column) and the corresponding potentials V((f>) (right column) for r ~ 0.18, n ~ 0.6, 
dnd/ log A; ~ —0.02. This parameter region is observationally disfavored, but shows the complicated behavior possible for 
solutions to the flow equations at higher order. 



recovered by a factor of 5 would allow one to start making 
meaningful reconstructions of the potential. 

V. DISCUSSION 

We have presented a novel method for reconstruct- 
ing the inflationary potential, based on the powerful 
flow equation approach to inflationary dynamics. This 
method, which we dub Monte Carlo reconstruction, dif- 
fers from other prescriptions for reconstruction in that it 
is stochastic in nature and involves extremely weak prior 
assumptions about the form of the potential. A stochas- 
tic approach to reconstruction has distinct advantages 
compared with, for example, methods which expand the 
potential in a Taylor series and fit the coefficients of the 
Taylor expansion to observable parameters. In particu- 
lar, the method can be applied equally well to both poor 
data sets and to high-quality data: parameters such as 
the tensor spectral index do not enter the reconstruction 
in any direct way, making the method very simple and 
robust. Also, Monte Carlo reconstruction naturally pro- 
duces shapes for the potential outside the region directly 
constrained by observation. This follows from our prin- 



cipal assumption, single field inflation, and the attractor 
behavior of the inflationary dynamics which tends to en- 
sure that the reconstructed potentials overlap outside the 
region in which we have direct observational input. 

We want to make clear what Monte Carlo reconstruc- 
tion is not: it is a stochastic method, but not a statistical 
one. We do not have a metric on the space of initial con- 
ditions. Consequently, we cannot use the "density" of 
models in any particular parameter space to infer the 
relative likelihood of one parameter region over another. 
Plots of models in either the space of possible potentials 
or the space of exponent m (from <fi m , as in Fig. 5) are 
properly interpreted as exclusion plots, indicating which 
regions are either consistent or inconsistent with the data. 
However, we cannot determine the relative likelihood of 
different initial points in slow roll space without an un- 
derstanding of inflationary initial conditions. Many other 
reconstruction attempts truncate the slow roll expansion 
at second or third order, but the approach here can be 
extended to arbitrary order in slow roll, and we have 
checked that our results do not depend on the specific 
level at which we truncate slow roll. However, any slow 
roll ansatz effectively rules out potentials with small fea- 
tures where the higher order derivatives of the potential 
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FIG. 4: The upper panel shows 100 reconstructed potentials, assuming r — 0.02 ±0.01, n = 0.95 ±0.01, dn/dlogk = 0.0 ±0.01 
(the errors bars expected from the Planck mission). This choice implies that the tensor modes are unambiguously detected, 
and leads to a tight constraint on the normalization of the potential. The field <f> is defined such that the observational 
parameters are calculated at (f> — 0. The lower plot shows 100 reconstructed potentials, for r — 0.0 ± 0.01, n = 0.93 ± 0.01, 
dn/dlog k = 0.0±0.01. In this case, tensor modes are not resolved, and the normalization of the potential is poorly constrained. 
The one anomalous potential in the top figure corresponds to a potential that, by chance, has comparatively large slow roll 
parameters, but which cancel in just the right way to produce the specified cosmological spectrum. 



and H((f>) are large p3| . If the feature is located in the 
region of the potential that corresponds to the primordial 
spectrum, strict limits can be placed on the size and slope 
of the feature [Q. However, if the feature is outside this 
region then we will obviously not be able to reconstruct 
it. 

The most straightforward application of Monte Carlo 
reconstruction is to simply generate an ensemble of po- 
tentials consistent with some observational constraint. 
This simple and robust procedure can be applied to data 
sets (for example the forthcoming data from the MAP 
satellite) for which direct functional reconstruction of the 
potential will likely be impossible. While this method 
suffers from the same fundamental limitations as any 
other method, it will allow us to at least answer qual- 
itative questions about the form of the inflationary po- 
tential: for example, is the potential convex or concave? 
Higher quality data will allow for more quantitative con- 
straints on the potential. Detection of a nonzero tensor 
component will give information about the normaliza- 



tion of the potential, to within an order of magnitude or 
two with the accuracy projected for Planck. Even ab- 
sent a detection of tensor modes, it may be possible to 
reach conclusions about the shape of the potential, if not 
its normalization. More quantitative information can be 
gleaned from more accurate data sets. For instance, an 
improvement in parameter resolution by a factor of five 
or so over Planck will make information about the ex- 
ponent <fi m of the potential available in at least a rough 
sense. Other possible applications would be plotting the 
results of the Monte Carlo reconstruction in the space 
of derivatives of the potential, V, V, V", and so on, for 
comparison with Refs. p5|, |2(|. 

Monte Carlo reconstruction is easily generalizable. For 
example, one need not use the slow roll approximation to 
calculate the power spectrum associated with a particu- 
lar choice of "slow roll" parameters generated by the flow 
equations. (Note that the slow roll expansion is com- 
pletely distinct from the slow roll approximation. The 
solutions we generate for the background evolution are 
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FIG. 5: These two histograms show the values of the power 
m (horizontal axes) obtained by fitting Eq. (^) to 100 po- 
tentials generated by the Monte Carlo reconstruction algo- 
rithm, where the measured spectra are assumed to have the 
central values predicted by A</> 4 inflation. The vertical axes 
indicate the number of models in a particular bin in m. In 
the top panel we assume that the error bars on the spec- 
tral parameters are equal to those expected from Planck, 
Sr ~ 0.01, Sn ~ 0.01 and Sdn/dlogk ~ 0.01, whereas the 
bottom panel corresponds to Sr ~ 0.002, Sn ~ 0.002 and 
Sdn/dlogk ~ 0.005. In the former case, the functional form 
of the potential cannot be meaningfully recovered, but in the 
lower case the results are consistent with V oc cj> 4 . In both 
plots we have dropped a few cases for which the least squares 
solver did not converge. 



exact.) For calculating the fluctuation spectrum associ- 
ated with a particular path in the parameter space, one 
could equally well apply the method of uniform approx- 
imations introduced by Habib, et al. |27). It is also 
straightforward to solve for the perturbation spectrum 
by numerically solving the exact equation for quantum 
modes in the inflationary spacetime: as in the calcula- 
tion of the potential itself, all the information required 
to do so is contained in the solution to the flow equations. 
Thus the method is not only exact in principle, but can be 
made so in practice if necessary. In addition, the same 
techniques could be applied to a set of flow equations 
based around an expansion other than the Hubble slow 
roll expansion of Liddle et al., such as the expansion used 
in Ref . j28) . Doing so would be useful to investigate the 
dependence of the reconstructed potentials on the details 
of the truncation scheme for the flow equations. We ex- 
pect this dependence to be small. 
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